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We have measured and analyzed the hyperfine structure of the only stable fermionic isotope of 
atomic erbium as well as determined its isotope shift relative to the four most-abundant bosonic 
isotopes. Our work focuses on the J ^ J + 1 laser cooling transitions from the [Xe]4f^'^6s'^(''Hei) 
ground state to two levels of the excited [Xe]4f^^6s6p configuration, which are of major interest for 
experiments on quantum-degenerate dipolar Fermi gases. From a fit to the observed spectra we 
find that the magnetic dipole and electric quadrupole hyperfine constants for the excited state of 
the strong optical transition at 401 nm are Ae/h = -100.1(3) MHz and Be/h = -3079(30) MHz, 
respectively. For the narrow transition at 583 nm a similar fit leads to hyperfine constants that 
are in excellent agreement with previous measurements. We have also determined the hyperfine 
constants using relativistic configuration-interaction ab-initio calculations. The agreement between 
the ab-mitio and fitted data for the ground state is better than 0.1%, while for the two excited states 
the agreement is 1% and 11% for the Ae and Be constants, respectively. 

PACS numbers: 32.10.Fn, 32.30.Jc, 37.10.De 
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I. INTRODUCTION 

The field of ultracold quantum gases has historically 
heavily relied on alkali-metal atoms. Only recently, the 
use of non-alkali-metal atoms has gained attention as a 
means to explore fascinating quantum phases of matter 
that are not accessible with alkali-metal species. Species 
with multiple unpaired valence electrons have rich atomic 
energy spectra and exhibit various types of coupling be- 
tween the electronic angular momentum J and the nu- 
clear spin / of the atom. For instance, fermionic alkaline- 
earth atoms have J = and 1^0 and the electronic and 
nuclear angular momenta decouple. This decoupling is at 
the center of proposals for efficient quantum simulation 
[1-3] and quantum magnetism [4-7]. Recently, degener- 
ate Bose and Fermi gases of Ca [8], Sr [9-11], and the 
alkaline-earth-like Yb atoms [12, 13] have been realized. 

Lanthanides with submerged 4f-shell electrons are a 
novel class of atoms that attract attention in the field of 
ultracold quantum physics. Lanthanide atoms can have 
an exceptionally large electronic angular momentum J 
resulting from the alignment of the angular momenta 
of the submerged electrons. Consequently, these species 
can have strong magnetic moments as large as 10 ^b, 
where /.is is the Bohr magneton. The mutual interac- 
tion is dominated by long-range magnetic dipole-dipole 
forces. Their dipolar character can be one hundred times 
larger than that for alkali-metal atoms. This key prop- 
erty makes lanthanides prime candidates for the study 
of atomic dipolar physics [14-16]. Two magnetic lan- 
thanide species, Dy [17, 18] and Er [19], have been re- 
cently brought to quantum degeneracy while others are 
under investigation [20-22]. 

The success of quantum-degenerate-gas experiments 



relies on a precise understanding of the atomic properties, 
such as energy levels, hyperfine structures, and atomic 
polarizabilities. However, for unconventional atomic 
species, such as lanthanides, the available knowledge is in 
many instances insufficient for laser cooling and trapping 
purposes. Therefore, dedicated experiments need to be 
conducted en route to quantum degeneracy [21, 23-25]. 

In this paper, we present a combined experimental 
and theoretical investigation of the hyperfine structure 
of the only stable fermionic erbium isotope, ^^^Er. In 
particular, we obtain the magnetic dipole. A, and elec- 
tric quadrupole, B, hyperfine structure constants for the 
ground and two electronically excited states of -^^'^Er, 
which are relevant for laser cooling experiments [23, 26]. 
The two electronic excited states investigated are the one 
at a wavelength of 582.67 nm (corresponding to photon 
energy E/{hc) = 17157.307 cm"!) and one at 400.796 nm 
{E/{hc) = 24943.272 cm~i) from the ground state [27]. 
Here h is Planck's constant and c is the speed of light. 
In addition to the study of the hyperfine constants, we 
also obtained the isotope shift of ^^^Er relative to the 
most-abundant bosonic isotopes. Our work provides im- 
portant information for future experiments on quantum- 
degenerate Fermi gases of strongly dipolar Er atoms. 

In a previous work, we used the optical transitions at 
about 401 nm and 583 nm for Zeeman slowing (ZS) and 
magneto-optical trapping (MOT) apphcations [26]. We 
demonstrated efficient laser cooling for five Er isotopes, 
including the fermionic one. However, the realization of a 
MOT of fermionic Er isotope was challenging since only 
the hyperfine structure of the ground and the 583 nm- 
excited state were known [28, 29], while the one of the 
state at 401 nm was unknown prior to this work. To oper- 
ate the Zeeman slower and the transversal cooling stage 



we had in fact to proceed empirically and try different 
locking points for the light at 401 nm before being able 
to produce a MOT of Fermions. 

Figure 1 shows the atomic level scheme of Er. The elec- 
tronic ground state belongs to the [Xe]4f^^6s^ configura- 
tion and has a large orbital angular momentum quantum 
number L — 5 (H-state) and a total electronic angular 
momentum quantum number J = 6. The excited states 
at 401 nm and 583 nm belong to the [Xe]4f^^6s6p config- 
uration and have singlet ^Pi and triplet '^Pi character 
for the outer two valence electrons, respectively. Both 
excited states have a total electron angular momentum 
J ~ 7. Erbium has six stable isotopes with natural abun- 
dance being 33.6 % for i^SEr, 26.8 % for ^^SEr, 23.0 % for 
is^Er , 14.9% for ^^OEr, 1.61% for i64Er, and 0.14% for 
i62gj. i67E]- ig tJie only stable fermonic isotope. The 
bosonic isotopes have zero nuclear spin (/ = 0) while the 
fermonic one has / = 7/2 and shows hyperfine structure. 
All three electronic states of ^^''Er have eight hyperfine 
levels ranging from F = J — 7/2 to _F = J -I- 7/2, where 

This paper is structured as follows. Section II de- 
scribes our experimental methods to investigate the hy- 
perfine structure of the relevant states. Section III re- 
ports on our least-squares fitting procedure to obtain the 
hyperfine constants and isotope shifts from the measured 
spectra. Section IV describes our ab-initio relativistic 
configuration-interaction calculations and compares the 
ab-initio hyperfine constants with the fitted values. Sec- 
tion V summarizes our conclusions. 
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FIG. 1: (color online). Energy levels of atomic Er up to 
E/{hc) — 25000 cm^^ for different electronic angular momen- 
tum quantum numbers J [27, 30]. States with odd (even) 
parity are indicated by black (red) horizontal lines. The two 
relevant laser cooling transitions at 401 nm and 583 nm are 
indicated by arrows. 



II. ATOMIC SPECTROSCOPY 

We measure the hyperfine structure of the ^^^Er iso- 
tope using modulation-transfer spectroscopy [31]. The 
spectroscopy is performed on an atomic Er vapor cre- 
ated with a hollow cathode discharge lamp (HCL). The 
HCL has the advantage to provide atomic vapors without 
the need of heating the source close to its melting point, 
which is very high in the case of Er (« 1529 °C). 

We use a commercially available HCL, which is filled 
with an argon gas at a fixed pressure of 4 mbar [32] . By 
applying a high voltage on the electrodes, the argon gas 
is ionized and accelerated into the center of the Er-coatcd 
cathode. When hitting the surface the kinetic energy of 
the Ar-ions is high enough to free neutral erbium atoms 
by sputtering processes [33]. We typically operate the 
HCL with a voltage of 110 V, giving a discharge current 
of 9.2 mA. 

We perform a Doppler-free modulation-transfer spec- 
troscopy in the HCL [34, 35] . The laser beam is split into 
a pump and a probe beam, as shown in Fig. 2(a). The 
pump light is modulated with an electro-optical modu- 
lator (EOM) driven by a local oscillator (LO) at a fre- 
quency of 14.2 MHz and with a power of 23dBm [36]. 
A four-wave mixing process transfers the sidebands from 
the modulated pump beam onto the counter-propagating 
probe beam [37]. We acquire the spectroscopy signal by 
mixing the LO signal with the probe beam signal de- 
tected by a photodiode (PD). By setting the LO and the 
signal either in-phase or shifted by tt, one can obtain a 
signal proportional to the dispersion or to the absorption 
of the atomic sample, respectively. In our setup we use 
the dispersion signal. 

In a first set of experiments, we measure the hyper- 
fine structure and the isotope shift of the excited state 
at 401 nm. A frequency-doubled diode laser is used for 
the spectroscopy. Figure 2(b) shows the dispersive spec- 
troscopy signal for this transition. The signal is averaged 
over 16 scans with a scanning speed of 2.4GHz/s. 

Our measurement reveals the full hyperfine structure 
for the fermionic ^^^Er. The discussion and the as- 
signment of the observed spectral features are given in 
Sec. HI. In addition, we determine the isotope shifts for 
the bosonic isotopes relative to ^^^Er. We measure a shift 
of -1681(14) MHz for i™Er, -840(14) MHz for ^^SEr, 
and +849(17) MHz for ^^^Er, which is in good agreement 
with Ref . [38] . The linewidths are extracted by fitting the 
derivative of a Lorentzian curve to the data. This gives an 
averaged value of 27r x 88(8) MHz, corresponding to about 
three times the natural linewidth of 2tt x 29.7(5) MHz 
[39] . This broadening of the transition can be explained 
as a combined effect of collisional and power broadening. 
For a number density of about 10^^ cm""^ and an argon 
background pressure of 4 mbar, we calculate a collisional 
broadening of 27r x 8.2 MHz. Considering a total intensity 
of the pump and probe beams of ^ = 250 mW/cm^, we 



estimate a power broadening of a factor of vl-i- '^ / >>^o 
2.3 with ^0 = 60.3mW/cm^ being the saturation inten- 
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FIG. 3: (color online). Spectroscopy signal and hyperfine as- 
signment of the 401-nm transition of the fermionic ^"^^Er. The 
black solid line is the recorded line. The red line is a simulated 
lineshape obtained using a nonlinear fit to the line positions 
and a linewidth of 7/(27r) — 90 MHz. The simulated lineshape 
is a sum of the first derivative of several Lorentzians, one for 
each hyperfine transition, whose relative strength is given by a 
theoretical estimate of the linestrength. We scaled the overall 
size of the simulated lineshape to fit to the experiment. The 
assignment of the P-branch transitions {Fg — s- -Fe — Fg + 1) 
is shown by vertical lines and pairs {Fg,Fe). The contribu- 
tion of the bosonic lines is not included in the theoretical 
simulation. The hyperfine coefficients of the excited state are 
Ae/h = -100.1 MHz and B^/h = -3079 MHz. 



FIG. 2: (color online). Modulation transfer spectroscopy of 
Er for the 401-nm and 583-nm transitions, (a) Laser setup 
for spectroscopy on a hollow cathode discharge lamp (HCL); 
see text. The pump (probe) light has a power of 3.3 mW 
(0.6 mW) for the 401-nm transition and 20 mW (ImW) for 
the 583-nm transition, (b) and (c) show the obtained spec- 
troscopy signals for the 401-nm and the 583-nm transitions of 
different isotopes. Signals related to the hyperfine structure 
of ^®^Er are indicated by arrows. The amplitudes of the ob- 
served signals are approximately proportional to the natural 
isotope abundances. 



sity. Combining the two contributions, we estimate a 
broadened linewidth of 2tt x 81 MHz, which is in agree- 
ment with the observed value. 

In a second set of measurements, we focus on the hy- 
perfine structure of the excited state at 583 nm, which has 
a linewidth of 27r x 186 kHz [40]. The spectroscopy is per- 
formed with a dye laser, which is frequency-stabilized to 
an ultra-low expansion cavity within 30 kHz [26] . We use 
a spectroscopy setup similar to the one described above 
for the 401-nm transition. Figure 2(c) shows the corre- 
sponding spectroscopy signal. 

Despite the narrow-line nature of the transition, we 
could observe five features related to the hyperfine struc- 
ture of the fermionic isotope and three features for the 
bosonic ones. The discussion of the hyperfine struc- 
ture is given in Sec. III. We measure an isotope shift of 



-975(15) MHz for ^^^F^t and -1966(14) MHz for i™Er 
relative to ^^^Er, respectively. These values are in good 
agreement with Ref. [29]. 

For this transition, we extract an averaged value for 
the linewidth of 2tt x 23(5) MHz, corresponding to about 
120 times the natural linewidth. This large broaden- 
ing can again be explained in term of coUisional and 
power broadening. Considering the saturation inten- 
sity of J^o = 0.13mW/cm^ and our total intensity of 
J^ := 1.3 x 10'^ mW/cm^, we calculate a power broadening 
of a factor of 100. Adding the effect of coUisional broad- 
ening, we obtain an overall linewidth of 27r x 19.3 MHz, 
which is in agreement with the measured value. Be- 
cause of this large broadening, we could operate the 
modulation-transfer spectroscopy at the same LO fre- 
quency as the one used for the 401-nm transition. 



III. ANALYSIS OF HYPERFINE STRUCTURE 

In this section we describe our fitting procedure to 
the observed spectra of the five most abundant Er iso- 
topes and we present the resulting hyperfine-structure 
constants A,, and B^ for ^^^Er. The bosonic features are 
easily assigned as shown in Fig. 2. The remaining weaker 
features, which sometimes overlap with those from the 
bosonic isotopes, are due to ^^^Er. 

We start with the definition of the transition energies 




detuning with respect to Er line (GHz) 

FIG. 4: (color online). Spectroscopy signal and hyperfine as- 
signment of the 583-nm transition of the fermionic ^ Er. The 
solid black line is the experimental spectrum while the red 
line is a simulated line-shape using a linewidth of 7/(27r) — 
20 MHz. The P-branch transitions {Fg -^ F^ = Fg + 1) are as- 
signed by pairs {Fg, Fe). Three P-branch resonances and sev- 
eral Q-branch [Fg — >■ Fe = Fg) resonances are predicted to lie 
outside of the measurement range. The simulated line-shape 
is a sum of the first derivative of several Lorentzians, one for 
each hyperfine transition, whose relative strength is given by a 
theoretical estimate of the line strength. We scaled the overall 
size of the simulated lineshape to fit to the experiment. The 
contribution of the bosonic lines is not included in the the- 
oretical simulation. The hyperfine coefficients of the excited 
state are A^/h = -f 72. 7 MHz and Be/h = -4457.2MHz. 



between the ground and an excited state of '^^^Er includ- 
ing hyperfine interactions [41] 

hhJF.Fg = Ai67 + tlhJiQfi + Ee{Fe, Je, I) ~ Eg{Fg, Jg, I) , 

(1) 

where Aig? is the ^^^Er isotope shift relative to the tran- 
sition energy fiwige of the bosonic ^^^Er atom, the most 
abundant isotope, and Ec,{Fe, Jg, /) and Eg{Fg, Jg, I) are 
hyperfine energies of the excited and ground state, re- 
spectively. The quantum numbers Fi and Ji with i = e 
or g are the total atomic and electronic angular momen- 
tum of the excited and ground state, respectively, and 



Ei{Fi, Ji,I) - 



3a(a + l)-4/(/+l)J,(J, + l) 
2/(2/ - 1)2 Ji(2J,-l) 



(2) 



where d = Fi{F^ + 1) - J^{J^ + 1) - /(/+ 1). Finally, the 
transition energies A^ -I- Swige define the isotope shift for 
bosonic Er isotopes with atomic number A. 

In addition to the resonance positions, we can calculate 
the line shape of the fermionic spectral features S{uj) by 
noting that the signal is well approximated by 



as a function of laser frequency w, where the sum is over 
all {Fg,Ff,) hyperfine lines, and C(ll),j) is a Lorentzian 
centered around zero with linewidth 7 [42]. Conse- 
quently, for an isolated line the resonance occurs when 
the signal is zero. The fiuorescence line strength Qf^,f 
is 

Qf..F, = Y. W9l)P9Mg\dlg\{JeI)F,Me)? 
M^Mgq 

\{JM\Je)\', (4) 



^^^^^45 5)) K'^^i""^^" 



S{uj) ex - ^ Q 



^-^'d:^(" 



^F,F„,l) 



(3) 



where the Mi are magnetic quantum numbers, dig is the 
electric dipole-moment operator, and we have assumed 
equal population for all hyperfine states F^Me of the 
electronic excited state. Finally, F = 2F + 1, ('■') is 
a six-j symbol, and (Jg||d|| Je) is a reduced dipole matrix 
element independent of Fg and F^ . 

We use a nonlinear least-squares fit to the experimen- 
tal spectra to determine the hyperfine constants and iso- 
tope shift Ai67 of the excited states. The fit is based on 
six resolved hyperfine features for the 401-nm line and 
five resolved features for the narrow 583-nm line. In our 
analysis, we hold the hyperfine constants for the ground 
state to the literature values of Ag/h = -120.487(1) MHz 
and Bg/h = -4552.984(10) MHz [28], which have sig- 
nificantly lower uncertainties than those for the excited 
states. 

Figure 3 and 4 are the results of our fit for the 
401-nm and 583-nm line, respectively. We observe re- 
markable agreements between the simulated and exper- 
imental spectra. For the excited 401-nm level, we ex- 
tract the best value for the hyperfine coefficients to be 
Ae/h = -100.1(3) MHz and Be/h = -3079(30) MHz. 
Using these coefficients and those for the ground state, 
we obtain resonance positions that agree to better than 
11 MHz with the experimental values. For the ex- 
cited 583-nm level, we fit the line shape of the spec- 
tral features while the resonance positions arc calculated 
by using the hyperfine constants of the excited states, 
Ae/h = -172.7MHz and Be/h = -4457.2MHz, from 
Ref. [43]. We note that the additional structure in the 
experimental data, which is not fitting to the theoretical 
curve, originates from a slightly misadjusted phase of the 
local oscillator in the spectroscopy setup. Table I com- 
pares the theoretical and experimental hyperfine energies 
Ai67 + Ee{Fe, JeJ) + Ee{Fg, Jg,I) for thc 583-nm and 
401-nm transitions in ^^^Er and lists the corresponding 
quantum numbers of Fg and Fe- Table II and Fig. 5 show 
the resulting isotope shifts A^ as a function of the mass 
number A relative to the energy of the ^^^Er isotope. 



IV. AB INITIO HYPERFINE CONSTANTS 

In conjunction with the experimental measurements 
and fits, we have performed extensive ab initio electronic 
structure calculations of the magnetic dipole A and elec- 
tric quadrupole B hyperfine constants. They describe 
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FIG. 5: (color online). Isotope shifts for tire 583-nm and 
401-nni lines of the isotopes ^^*Er up to '^^"Er as a function 
of mass number where the transition energy for the bosonic 
isotope ^ Er is taken as energy reference. The isotope shift 
of the bosonic isotopes falls on a single straight line, with the 
isotope shift of the center of gravity of the fermionic ^''^Er 
isotope, green cross and green square, is slightly displaced 
from this linear dependence. 



constants have been performed using a relativistic mul- 
ticonfiguration Dirac-Fock (MCDF) method [44]. In this 
method we perform an all-electron calculation of the 
wave function leading to an accurate description of the 
electron-spin density near the nucleus. The eigenfunc- 
tions are superpositions of non-orthogonal many-electron 
determinants of one-electron Dirac-Fock functions for the 
core and valence orbitals and Sturm functions for virtual 
orbitals. Both types of one-electron orbitals are opti- 
mized for either the Af^'^6s'^ ground or 4/^^6s6p excited 
state reference configurations. 

The hyperfine splittings of atomic levels are due to in- 
teractions between electrons and nuclear multipole mo- 
ments. In the configuration interaction picture and using 
atomic units the A and B constants are given by 



gifJ-N 
Mj 



A{J) 



B{J) - 7^ 



(v,,JA/,|^llii|lk|^,jM,),(5) 




2J(2J-1)(2J-H) 

(2J + 2)(2J + 3) 



(6) 



the coupling of the nuclear spin / to the total electron 
angular momentum J, due to the magnetic dipole and 
electric quadrupole interaction, respectively. The latter 
originates from the electric field gradient created by the 
electrons at the nuclear location. We were interested to 
reproduce the known constants for the Er ground state 
as well as those of the excited state at the 583-nm line 
obtained by Ref. [28, 43]. We can then confirm our mea- 
surement of the unknown constants of the excited level 
at the 401-nm line. 

The ah initio calculations of the hyperfine structure 



TABLE I: The observed and calculated hyperfine energies 
Ai67 + Ee{Fe, Je , /) + Ee{Fg, Jg , I) for the 583-nm and 401- 
nm lines in ^®^Er. The theoretical values are based on 
the hyperfine coefficients Ae/h = —172.7 MHz and B^/h = 
-4457.2 MHz [43] for the 583-nm line and our values Ae/h = 
-100.1 MHz and Be/h = -3079 MHz for the 401-nm line. 
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where the sum i is over all electrons with positions fi 
with respect to the nucleus, Ye„i{f) are spherical har- 
monics, and di is the Dirac matrix for electron i. Fur- 
thermore, gi is the nuclear g-factor, /z^v is the nu- 
clear magneton in atomic units, and Q is the nuclear 
quadrupole moment. The relativistic electronic eigen- 
functions |\1/, JMj) = ^pCp\(j)p,JMj), obtained from 
the configuration-interaction calculations with relativis- 
tic determinants \4ip,JMj) and CI coefficients cp, have 
total angular momentum J and projection Mj. 

When an atom has open or unfilled electron shells, it 
leads to an unbalanced electron-spin density near the lo- 
cation of the nucleus. As hyperfine constants are propor- 
tional to the difference in electron spin densities this leads 
to non-zero A and B coefficients. To account for this ef- 
fect we use a model, where the single-electron orbitals 



TABLE II: The observed isotope shift of the Er isotopes for 
the 583-nm and 401-nm lines. The transition energy for the 
bosonic isotope ^^^Er is taken as energy reference. The iso- 
tope shift of the center of gravity of fermionic ^^^Er was ob- 
tained from fitting. 
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Energy 
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583 


nm 


401 


nm 


170 


-1966(14) 




-1681(14) 




168 


-975(15) 




-840(14) 




167 




-337(11) 




-297(6) 


166 
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- 




+849(17) 





TABLE III: A summary of the relevant hyperfine A and B 
constants for ground and excited states (e.s.) usable for laser- 
cooling of ^''^Er. 



state 


J 


A/h (MHz) 


B/h (MHz) 


Rcf. 


ground state 
ab initio 


6 


-120.487(1) 
-120.42 


-4552.984(10) 
-4554 


[28] 
this work 


583-nm e.s. 
ab initio 


7 


-172.70(7) 
-174 


-4457.2(29) 
-4057 


[43] 
this work 


401-nm e.s. 
ab initio 


7 


-100.1(3) 
-100 


-3079(30) 
-3424 


this work 
this work 



differ for each spin direction or more precisely for each 
spinor of the Dirac-Fock equation. Alternatively, this im- 
plies different exchange potentials for electrons with spin 
up or down. 

We use three restricted active spaces (RAS) to clas- 
sify the electron Dirac-Fock and Sturm orbitals, ensur- 
ing an efhcient and compact CI expansion that, never- 
theless, remains accurate. The first group of orbitals, 
RASl, contains the occupied spinors of the relevant ref- 
erence configuration. We have studied convergence of 
the hyperfine structure constants as the active set of 
orbitals was systematically increased. For our most- 
precise Er atom calculation RASl includes the occupied 



Is ,2s , 2p^^2J '^P3/2i 



.., and 4/|,2 shell electrons. We al- 



low up to one electron to be excited out of RASl into the 
two other active spaces. The second group, RAS2, con- 
tains the open-shell 4/^,2, ^^1/2^ ^^^ ^Pi/2, 6^3/2 spinors, 
while the third group, RAS3, contains spinors that are 
unoccupied in the reference configuration. These latter 
virtual Sturm orbitals are the high-lying s-wave spinors 
from 7s up to 13s, p-wave spinors from 7p up to lip, and 
the 5(i-spinor. For both RAS2 and RAS3 we allow up to 
two electrons to enter or leave. 

With this basis our finite-nuclear-size and finite- 
nuclear mass corrected ab initio values of the A and B 
hyperfine constants are —120.42 MHz and —4554 MHz 
for the ground state level, -174 MHz and -4057 MHz 
for the excited level at 583 nm, and —100 MHz and 
—3424 MHz for the excited level at 401 nm, respec- 
tively. Consequently, the ab initio A constants agree 
with experimentally-determined values to better than 
1%, whereas the B constants differ by up to 11% for the 
two excited states. For the ground state the agreement 
for the B constant is also better than 1%. 



levels within the excited 4/^^6s6p configuration. A 
least-squares algorithm applied to the experimentally- 
measured hyperfine-structure energies gives accurate val- 
ues for the two constants as well as values for the isotope 
shift of five isotopes. The ab initio calculation is based 
on a multi-conflguartion Dirac-Fock method where we al- 
low no more than two electrons to be excited from and 
between the active spaces. The method has no further 
adjustable parameters. 



Our results are summarized in Table III and Fig. 6. We 
find that the ab initio A coefficients for all three states 
and the B coefficient for the ground state agree to bet- 
ter than 1% with the experimental values, which is in 
a surprisingly good agreement considering the complex 
electron-shell structure of the Er atom. We note that 
the ab initio electric quadrupole constants B for the two 
excited states exhibit a larger deviation from the exper- 
imental values. This might be a consequence of missing 
key configurations: the excited states have three open 
shells, 4/^^, 6s, and 6p, from which more than two elec- 
trons might need to be excited. In addition, Sternheimer 
shielding (e.g., distortions in the electron shells by the 
nuclear quadrupole moment), which is not considered in 
our MCDF theory, might cause significant corrections. 



4f 26s6p (ip,) 




^Er 



4f26s6p (3P^) 



25.8 MHz 

■ 135MHz 

■ 463 MHz 

■ 902 MHz 

■ 1.50 GHz 

■ 2.27 GHz 
■3.25 GHz 



69.1 MHz 
234 MHz 
44.5MHz 
354 MHz 
468 MHz 
1.65GHz 
2.69GHz 



V. CONCLUSIONS 

We have used laser modulation-transfer spectroscopy 
on atomic Er as well as performed ab initio electronic 
structure calculations of Er to obtain the magnetic dipole 
and electric quadrupole constants for the only stable 
fermionic isotope, ^^^Er. We focusscd on transitions 
from the Af^^ds^ (J — 6) ground state to two J = 7 



FIG. 6: (color online). Hyperfine levels of the ground (g.s.) 
and two excited states of ^^^Er with particular interest for 
laser-cooling. The level splitting was calculated using A and 
B constants given in Table III for the respective transitions. 
The arrows depict two laser-cooling transitions. The transi- 
tion at 401 nm used for Zeeman slowing is shown in blue and 
the transition used for magneto-optical trapping at 583 nm is 
shown in yellow. 
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